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ABSTRACT 

In searches for low-mass companions to late-type stars, correlation between radial velocity variations 
and line bisector slope changes indicates contamination by large starspots. Two young stars demonstrate 
that this test is not sufficient to rule out starspots as a cause of radial velocity variations. As part of our 
survey for substellar companions to T Tauri stars, we identified the ^2 Myr old planet host candidates 
DN Tau and V836 Tau. In both cases, visible light radial velocity modulation appears periodic and is 
uncorrelated with line bisector span variations, suggesting close companions of several M.j up in these 
systems. However, high-resolution, infrared spectroscopy shows that starspots cause the radial velocity 
variations. We also report unambiguous results for V827 Tau, identified as a spotted star on the basis 
of both visible light and infrared spectroscopy. Our results suggest that infrared follow up observations 
are critical for determining the source of radial velocity modulation in young, spotted stars. 

Subject headings: planetary systems: formation - stars: individual (DN Tau, V836 Tau, V827 Tau) - 
stars: spots techniques: radial velocities 
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1. INTRODUCTION 

Extrasolar planets are common; over 300 systems have 
been discovered (e.g., ex oplanet.eu) . Recent studies have 
target ed higher m ass (|Sato et all l2007t Johnson et al.l 
l2007af) . lower mass dButler et al.ll2006bl: lEndl et al.ll200lT 
and y ounger objects (jPaulson et al . 2006; Setiawan et al.l 
120070 . Identifying young planets is important to define 
the time scale for planet formation and thus distinguish 
the possible formation process(es). 

Young stars still surrounded by the circumstellar ma- 
terial which forms planets are typically located at dis- 
tances of >100 pc and are thus inherently faint and of- 
ten o bscured. They also manifest strong magnetic activity 
(e.g.. lJohns-Krul]||2007D and are highly spotted. Numer- 
ous, large spots complicate detection of extrasolar planets 
throu gh radial velocity (RV) monitoring (|Saar fc Donahue] 
Il997f) because a spot that is partially visible at all times 
on th e surface of an inclined star mimics RV modulation 
(e g.jBouvier et a.1.1 120071 : iHuerta et al.ll2008D. 

iPaulson et alj (|2006f) studied 12-300 Myr old nearby 
stars and found no evide nce for planets with masses 
>l-2 M Jup at the 3 a level. ISetiawan et all (|2007f) identi- 
fied a minimum mass 6.1 Mj up planet in a 852 day period 
orbit around t he 100 Myr old G1-G 1.5 V star HD 70573. 
More recently. ISetiawan et alj (|2008l ) reported a ~10 Mj up 
planet in a 3.5 6 day orbit around th e 10 Myr old star TW 
Hya, although lHuelamo et al.l ([2008) identify this result as 
attributable to spots. 

In this letter we present our observations of the young 
stars DN Tau, V836 Tau, and V827 Tau. While the V827 
Tau visible light data clearly implicate spots as the cause 
of the apparent RV variability, corresponding data for DN 



Tau and V836 Tau suggest the presence of giant planets. 
Our infrared (IR) observations show that spots cause the 
RV variations seen in all three stars. In §2 we describe 
the observations, in §3 present our data analysis and the 
evidence for spots, and in §4 provide a brief discussion. 
We summarize in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Visible Light Spectroscopy 

Visible light spectra of DN Tau (M0), V836 Tau (K7), 
and V827 Tau (K7) were taken at the McDonald Oberva- 
tory 2.7 meter Harlan J. Smith telescope, between Novem- 
ber, 2004 and January, 20 08, with the Coude echelle spec- 
trograph (jTull et al.lll995h . A 1.2" slit yielded R-60,000. 
Integration times were ^1800 s; average seeing was ~2". 
ThAr exposures taken immediately before and after each 
spectrum provided wavelength calibration; typical RMS 
values for the di spersion solution prec i sion were ~4 m s -1 . 
RV standards (iNidever et all l2002t iButler et all 119961: 
ICumming et alj fl999) were observed on every night of ev- 
ery run; their overall RMS scatter is 140 m s -1 . We ob- 
tained 43 spectra of DN Tau, 21 of V836 Tau, and 20 
of V827 Tau and applie d standard IRA F r eduction rou- 
tines. Details are given in lHuertal (|2007f ) and lHuerta et al.1 

2.2. Infrared Spectroscopy 



We observed DN Tau, V836 Tau, V827 Tau, and the 
RV standards HD 65277 (K5) and G J 281 (MO) on 
UT 2008 February 13-2 with CSHELL (jTokunaga et all 
[19901 lG reene et al]|1993[ ). the high- resolution, IR spectro- 
graph at the NASA IRTF 3-m telescope. The seeing was 

1 This paper includes data taken at The McDonald Observatory of The University of Texas at Austin. 

2 Lowell Observatory, 1400 West Mars Hill Rd. Flagstaff, AZ 86001; lprato@lowell.edu 

3 Visiting Astronomer at the Infrared Telescope Facility, which is operated by the University of Hawaii under cooperative agreement NCC 
5-538 with the National Aeronautics and Space Administration, Office of Space Science, Planetary Astronomy Program. 

4 American Astronomical Society, 2000 Florida Avenue Suite 400, Washington, DC 20009; marcos.hucrta@aas.org 

5 Department of Physics and Astronomy, Rice University, MS-108, 6100 Main Street, Houston, TX 77005; cmj@rice.edu, naved@rice.edu, 
hartigan@rice.edu 

6 Department of Astronomy, University of Texas, R. L. Moore Hall, Austin, TX 78712; dtj@astro.as.utexas.edu 

1 



2 



0.4—0.8". Each object was observed on 6—8 nights. The 
0.5" slit yielded R~46,000. We obtained data in 10" nod- 
ded pairs. Spectra were centered at 2.298 fim. (vacuum). 
Integration times were ~1 hour for the T Tauri stars and 
—8 minutes for the standards. The signal to noise ratio 
(SNR) was —70— 120. D ata were reduced as described in 
Uohns-Krull et alJ (1 19991) . 

3. ANALYSIS 
3.1. Radial Velocities: Visible Light 

Relative RVs were determined by cross correlating a 
high SNR, fiducial spectrum against all other spectra for 
the same target. We used six orders spanning —5700 A to 
6800 A. Uncertainties were estimated from the standard 
deviation of the mean for the 6 orders, added in quadra- 
ture with the 140 m s -1 uncertai nty derived from the RV 
standards (52; lHuerta et al. 200gj). RVs were corrected for 
the Earth's barycentric motion. 

Optimum periods and uncertainties for p hasing the RV 
data were selected based on power spectra ()Huerta et al.l 
120081) . For DN Tau we found P=6.33 ±0.20 days and a 
false alarm probability (FAP) of <0.001, for V836 Tau, 
P=2.48 ±0.49 days and FAP=0.10, and for V827 Tau, 
P=3.76 ±0.06 days and FAP<0.001. We also checked 
for periodicity using the discrete Fourier transform plus 
CLEAN method of Roberts et al. (1987). The strongest 
power spectrum peaks for DN Tau and V836 Tau occur at 
the same periods. The CLEAN method recovered a best 
period of 3.61 d for V827 Tau, within —2 a of the above 
estimate. The phased RV data are shown in Figures 1—3. 

3.2. Line Bisector Spans 

The presence of a starspot will distort the line profile 
at the RV that corresponds to the stellar velocity at the 
location of the spot. This distortion is in proportion to 
the ratio of quiescent photosphere surface brightness and 
surface brightness within the spot, at the observing wave- 
length, an d the fraction of st ellar surface covered by the 
spot (e.g.. iQueloz et al1l2001[ ). Thus, asymmetries in the 
line profiles originating from spots will be present for all 
lines in the spectrum of a young star and are typically 
correlated with the RV measured from the same spec- 
trum. These asymmetries have become the standard crite- 
rion for rejecting starspots as the cause of false RV signals 
(e.g., IQueloz et alJl2001llBouvier et al.|[2007t iHuerta et al.l 
20081 ISetiawan et al.ll2008f) . For each of the six orders used 
to determine the RVs, we cross-correlated all absorption 
lines and measured the cross-correlation function (CCF) 
for that order. The average of these six CCFs was used 
to measure the bisector spans (lower panels of Figures 
1—3). The linear correlation coeffi cient and associated 
FAP (|Bevington fe Robinson! |1992| ) is listed in the cap- 
tions. As expected for a spotted star, a clear correlation 
between bisector span and RV is observed for V827 Tau. 
DN Tau and V836 Tau show no correlation, suggesting 
that the variability is not the result of spots. 

3.3. Infrared Radial Velocity Signatures 

The contr ast between a 4000 K photosphere and a 
3000 K spot (H ouvier fe Bertoudll989l) is greater in visible 
than in IR light because flux scales as a steeper function of 



temperatur e at wavelengths short er than the black body 
peak (e.g.. ICarpenter et al.l l200lh . Given the decreased 
spot to photosphere contrast in the near-IR, the ampli- 
tude of the RV modulation will be smaller. Conversely, if 
a planet drives the RV modulation, the amplitude should 
b e the same in visible and IR light. 

iBlake et ail (|2007l [2008) used near-IR observations to 
search for companions to low-mass objects, exploiting tel- 
luric absorption lines for high-precision RV measurements. 
Figure 4 shows an example of a GJ 281 K-band spec- 
trum and illustrates our similar approach. We created 
models by combining high resolution telluric absorption 
(jLivingston fe Wal lace 199l|) and cool stellar spectra (the 
sunspot atlas of Wallace & Livingston 1992), applying a 
range of velocity shifts relative to the telluric lines. Other 
free parameters are vsini, a Gaussian FWHM for the 
spectrometer line spread function, scale factors for line 
depths, and a first order continuum normalization func- 
tion. We employed the Marquardt method for no n-linear 
least squares fitting (|Bevington fc Robinson|[T992l ) of each 
model to an observed spectrum. The difference between 
the stellar and telluric velocities in the best fit model spec- 
trum yields the RV, which was then corrected for barycen- 
tric motion. 

Figures 1—3 show the IR-derived relative RVs. The data 
are phased to the periods given in §3.1. The RV standard 
deviation in the IR data for HD 65277 is 127 m s _1 and for 
GJ 281 is 98 m s" 1 . Internal errors, measured from the 
least squares fitting, are —40 m s _1 for both standards. 
For the young stars, internal errors were 100—300 m s _1 , 
depending on the SNR achieved. Random errors, which 
we assume add in quadrature with our internal errors to 
give the overall scatter in velocities, are 120 m s _1 for 
HD 65277 and 90 m s" 1 for GJ 281. We use 110 m s" 1 
as our final value for the random errors. The uncertain- 
ties shown for the IR data in Figures 1—3 represent the 
sum, in quadrature, of the individual internal error and 
the 110 m s _1 random error. Within our measurement 
precision, we are unable to detect any IR RV variability 
in DN Tau and V836 Tau. V827 Tau shows significant 
IR RV variations but at a reduced amplitude from those 
observed in visible light. 

4. DISCUSSION 

Figures 1 and 2 show the visible light RV modulation for 
DN Tau and V836 Tau, with full amplitudes of ~ 1500 and 
~ 2700 m s _1 , respectively. Within the 1 a uncertainties, 
all but one point in the IR RVs of DN Tau are consistent 
with zero; all six IR RVs of V836 Tau also show no vari- 
ation. These results indicate that no planets are present 
around DN Tau or V836 Tau with masses greater than a 
few M Jup at <0.5 AU or —10 M Jup at —1 AU, despite the 
absence of a cor relation betwee n the v isible light RVs and 
bisector spans. iMathieu et alJ (|1989l ) identify V836 Tau 
and V827 Tau as RV variables with peak-to-peak ampli- 
tudes of 7—8 km s _1 . Apparently the density and size of 
spots on V836 Tau vary; historical data may serve as an 
additional criterion for heavily spotted young stars. 

The primary conclusion from our visible light data is 
that the lack of a correlation between the line bisector 
span and the RV is not proof of a reflex motion compan- 
ion. In addition, the RV period can change significantly 



3 



with new data. Initial analysis of 20 visible light RVs 
for DN Tau, taken over 2.5 years, indicated convincing 
modulation with P=7.5 days and FAP=0.002. Although 
not markedly different from rotat i on period estimat es of 
6.0-6.4 days (|Bouvier et all 119931; IPercv et al.ll2006l ). x 2 
for P=7.5 days was more favorable than that for a sec- 
ondary peak at 6.3 days. Removing the fit to either the 
7.5 or 6.3 day period from the RV data and recalculating 
the power spectrum yielded no significant peak, suggesting 
that the data were best represented with only one period. 
With the addition of RVs measured in winter 2007-2008, 
the integrity of the 7.5 day period was diminished and the 
6.3 day period came to dominate the power spectrum. The 
rotation peri od of V836 Tau has been stable for decades 
at 6.76 days (|Grankin et al.ll2008l) . We find a different RV 
period, 2.48 days, and substantially reduced power in the 
RV modulation near 6.76 days. The RV period for V827 
Tau, 3.76 days, is indistinguishable f rom previously deter - 
mined values of the rotation period (jGrankin et a l. 2008). 

If the visible light RV modulation originates in spots, 
why don't the bisector spans correlate (Figures 1 and 2)1 
Figure 3 clearly shows correlation for V82 7 Tau; many 
other stars show a correlation as well (e. g., Queloz et all 



120011; iBouvier et alj 120071 : iHuerta et al.ll2008f ). We do not 

believe that this can be attributed to a stronger impact of 
spots on particular spectral lines because the same echelle 
order s were used to determ ine the bisector spans for all tar- 
gets. iDesort et al.l (|2007l ) show that when vsmi is smaller 
than the spectrometer resolution, RV and bisector varia- 
tions originating in spots can mimic the behavior expected 
from short period giant planets. The vsird of DN Tau and 
V836 Tau, ^10 km s _1 , are not much larger than our vis- 
ible light spectral resolution, 5 km s _1 , while V827 Tau 
and LkCa 19 (Huerta et al. 2008) have usini values of 
~20 km s _1 , suggesting qualitat ive agreement with the 
simulations of lDesort et alj (|2007[ ). 

5. SUMMARY 

We have measured the RVs of the ~2 Myr old T Tauri 
stars DN Tau, V836 Tau, and V827 Tau in visible and IR 



light; the variations we see in all three systems are likely 
the result of starspots. Furthermore: (1) Periodogram 
analysis can reveal the presence of a period but not nec- 
essarily a reliable value until RV measurements densely 
sample the full phase of the periodic signal. (2) The lack 
of correlation of line bisector spans with RVs for some stars 
with spots is an important problem to address in the search 
for planets around young stars. (3) High-resolution, IR 
spectroscopy is critical for the verification of young planet 
candidates; without the contrast in RV modulation ampli- 
tude between the visible light and IR data, it is unclear 
whether spots or a companion cause variability. (4) In the 
case of the spotted young star, V827 Tau (Figure 3), we 
observe a substantial reduction in RV amplitude between 
the visible light and IR data. This result can be exploited 
to improve our understanding of starspot temperature and 
filling factors. (5) T Tauri stars are virtually guaranteed 
to have spots at some level; if planets are also present, 
their detection will likely require disentanglement of the 
blended RV signals and will benefit greatly from measure- 
ments that span a broad range in wavelength. 
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Fig. 1. — Top: Relative RV versus phase for DN Tau. Visible light (open squares) and IR (filled circles) data have been phased to a priod 
of 6.33 ±0.20 days. Bottom: Line bisector spans as a function of RV for the visible light data. The linear correlation coefficient is r = —0.27 
and the associated FAP, f p , is 0.09. 
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Fig. 2. — Same as Figure 1 for V836 Tau, phased to a period of 2.48 ±0.49 days (r = 



-0.30, f p = 0.20). 
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Fig. 3. — Same as Figure 1 for V827 Tau, phased to a period of 3.76 ±0.06 days (r = 



-0.88, f p = 2.8 x 10~ 7 ). 
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Fig. 4. — Middle: UT 20 February 2008 CSHELL spectrum (solid line) of the RV standard GJ 281 and a model (dash-dot line) consisting 
of a stellar photosphere component (Top: from the NSO sunspot atlas) and a telluric component (Bottom: from the NOAO spectral atlas). 



